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I. INTRODUCTION 

One-dimensional directed polymers in a quenched random potential and equivalent problem of the solutions of 
the KPZ-equation Q 

describing the growth in time of an interface in the presence of noise have been the subject 
of intense investigations during the past two decades (see e.g. 0-0] )■ The model of directed polymers describes an 
elastic string directed along the r-axis within an interval [0, i]. Randomness enters the problem through a disorder 
potential V[4>(t),t], which competes against the elastic energy. The system is defined by the Hamiltonian 

ft 

H 



)>Vl=l dr{i[a r 0(r)] 2 + y[^(r),r]}; (1) 



where the disorder potential V[<^, t] is Gaussian distributed with a zero mean V(<p,r) — and the ^-correlations: 



V(<j), t)V(4/, t') = uS(t - T r )6{<j) - 0') (2) 

Here the parameter u describes the strength of the disorder. 

In what follows we consider the problem in which the polymer is fixed at the origin, <fi(0) = and it is free at r = t. 
In other words, for a given realization of the random potential V the partition function of the considered system is: 



Z = dx Z(x) = expl-/?^} (3) 

J — oo 

where 

Z(x) = / Xty( r ) Q- f3H ^ (4) 

J0(O)=O 

is the partition function of the system with the fixed boundary condition, <f>(t) = x and F is the total free energy. 
Besides the usual extensive part f$t (where /o is the linear free energy density) , the total free energy F of such system 
is known to contain the disorder dependent fluctuating contribution which in the limit of large t scales as t 1 ^ 3 (see 
e.g. [H3). 

In other words, in the limit of large t the total (random) free energy of the system can be represented 
as F — fat + at 1 / 3, f ', where c is a non-universal parameter, which depends on the temperature and the strength of 
disorder, and / is the random quantity which in the thermodynamic limit t — > oo is described by a non-trivial universal 
distribution function P(f). The trivial self-averaging contribution fot to the free energy can be eliminated from the 
further study by the simple redefinition of the partition function, Z = exp{— /3/oi} Z, so that Z = exp{— A/}, where 
A = /Sci 1 / 3 . Thus, to simplify notations the contribution fot will be just dropped out in the further calculations. 

For the problem with the zero boundary conditions, 0(0) = <p{t) = 0, the distribution function P(f) was proved to 
be described by the Gaussian Unitary Ensemble (GUE) Tracy- Widom distribution f§-fll"j|. O n the other hand, the free 
energy distribution function of the directed polymers with the free boundary conditions, eqs.([T])-(|4]), was shown to be 
given by the Gaussian Orthogonal Ensemble (GOE) Tracy- Widom distribution [l2|, [l3[ In the course of these proofs 
rather efficient Bethe ansatz replica technique has been developed [I3-[l3j . Here in terms of this technique we are going 
to study on the statistical properties of the transverse fluctuations of the directed polymers. The scaling properties of 
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the typical value of the endpoint deviations, 4>(t), et large times is well known: {<j)(t)) 2 oc i 4 / 3 (here (...) denotes the 
thermal average and (...) is the average over the disorder potential, eq.©) Qt3|- Much more interesting object is the 
probability distribution function P(x) for the rescaled quantity x = <f>(t)/t 2 ' 3 which is expected to become a universal 
function in the limit t — > oo. Recently this function has been derived in terms of the so called maximal point of the 
Airy 2 process minus a parabola jl4Ml6| . which is believed to descibe the scaling limit of the endpoint of the directed 
polymers in a random potential. The obtained explicit expression for P{x) turned out to be rather complicated and 
its analytic properties is not so easy to analyze although the asymptotic behavior of this function is already known: 
P(x -> oo) ~ exp{-|a;|7l2}[l5j]. 

In this work the explicit form of the distribution function of the directed polymer's endpoint fluctuations will be 
derived in terms of the Bethe ansatz replica technique. The distribution function we are going to consider is defined 
as follows: 



W(x) 



lim Prob[(j)(t)t 



-2/3 



> x\ 



dx P{x') 



(5) 



This function gives the probability that the rescaled value of the polymer's right endpoint <p{t)/t 2 ^ is bigger than a 
given value x. In this paper it will be shown that (see eas. (|79|) - (l84|) below) 



W(x) 



+ 00 f f + OG * \~1 

dfF 1 {-}) cLjcLv'U-bA * u , w (/,x) 

./ ./n \ / oioi' 



(6) 



where i r i(— /) = det[l — -B-/1 is the GOE Tracy- Widom distribution with the kernel B-f(ui,uj') 
/), (uj,uj' > 0) and 



Ai(w 



$u,u,> {f, x) 



+ 00 



dy 



\duj 



_d 



7 )*(w --/ + »; a;)*( w '-i/ + y; -x) 



where the function of two variables ^(uj;x) is the simple generalization of the Airy function: 



ty(oj;x) 



f dz rl 3 1 2 -i 

/ t: — exp-^ —z — tuz — —xz > 



(7) 



(8) 



Here the integration contour C in the complex plane starts at infinity with the argument — 7r/3 and ends up at infinity 
with the argument +7r/3. The above result looks quite similar to the one obtained in [14|, although at the moment I 
am not able to provide the proof that these results are indeed the same. 

The paper is organized as follows. In Section II we define the distribution function W(x) via the two-point free 
energy distribution function V x (fi, f%) which give the probability that the free energy of the polymer with the endpoint 
located above a position x is bigger than a given value /i, while the free energy of the polymer with the endpoint 
located below the position x is bigger than a given value ji . In Section III the function V x {J\^f%) is defined by mapping 
the considered problem to the one-dimensional A-particle system of quantum bosons with attractive ^-interactions. 
In Section IV the explicit expression for the probability function Vx(j\\fz) is obtained in terms of the Bethe ansatz 
replica technique. Finally, in Section V the result eqs.([5])-([5|) is derived. Conclusions and future perspectives are 
discussed in Section VI. 



II. THE ENDPOINT PROBABILITY DISTRIBUTION FUNCTION 



In terms of the partition function Z(x), eqs.@, the probability distribution function of the polymer's endpoint 
W(x), eq.©, can be defined as follows: 



in,; := Inn j / 



t— >oo 



lim { ) (9) 
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where 



Z { -\x) = f dx' Z(x') = exp{-A/ ( _)} (10) 

J —oo 



+00 



Z^\x) = / dx' Z(x') = exp{-A/ (+) } (11) 



where the parameter A oc t 1 / 3 and are the free energies of the polymers with the endpoint cf>(t) located corre- 
spondingly above and below a given position x. According to these definitions we find 

w< y v cxp{-A/ (+) } / °. fOT A-) < A+) 

W{x) = hm — — — — — — — — - = < (12) 

a^oo cxp{-A/ ( _)} + cxp{-A/ (+) } y 1; for /( ) > /(+) 

Let us introduce the joint probability density function V x [/(+) >/(-)] • By definition the quantity 
"Px [/(+)! /(-)] df(+)df(-) gives the probability that the free energy of the polymer with the endpoint located be- 
low x is equal to (within the interval d/(_)), while the free energy of the polymer with the endpoint located above 
x is equal to f/ + \ (within the interval <i/(+)). Thus, according to eq.flTT" 



f+00 p+00 

W(x) - / df (+) / [/(+)!/(-)] (13) 

1 -'/(+) 



Let us introduce one more joint probability distribution function: 

p + co f+00 

V x (fi,h) = Prob[/ (+) >/ i; / ( _) >/ a ] = / <#(+)/ 4f ( _)P B [/(+);/(_)] 



(14) 



This two-point free energy distribution function gives the probability that the free energy of the polymer with the 
endpoint located above the position x is bigger than a given value /1, while the free energy of the polymer with the 
endpoint located below the position x is bigger than a given value ji. According to this definition, 

V x [h\h\ = ^ r ^ r V x (f l ,f 2 ) (15) 
oh 072 

Substituting this relation into eq. ([T3")) we find 

W(x = / dh df 2 — — V x (h,f 2 ) (16) 

7-oo Jf 1 ofi df 2 

Integrating by parts over f 2 and taking into account that V x (fi, f 2 ) \, _ , =0 we get 

W(x) = - f + °° df 1 (^-V x (fij2)) (17) 

J-OO V OJl ' /2=/l+0 

Thus, to get the distribution function W{x) for the polymer's endpoint fluctuations we have to derive the two-point 
free energy distribution function V x (f\,f 2 ) first. Note that this function is different from the two-point free energy 
distribution function derived in [13] which describes joint statistics of the free energies of the directed polymers coming 
to two different endpoints 

III. MAPPING TO QUANTUM BOSONS 

According to the definition, eq. (fT4| . the probability distribution function V x (fi,f 2 ) can be defined as follows: 



00 00 



V x (h,f 2 ) = hm E ( li )L( ^ )fl exp(A£/i+Afl/ 2 ) [ZW(x)] L [Z(-)(x)] R (18) 
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Indeed, substituting here the definitions, eqs. (ITOl) - (fTT|) , we find: 



V x {h,f 2 ) 



/+oo r+oo 00 00 (_-\\L (_\\R 

df i+) / [/(+)? /(-)] E E S^Sr ex P[ AL (/i - /(+))] exp[Ai?(/ 2 - /(_))] 



L=0 i?=0 



lim 

A— ^00 



+00 



+00 



d/ (+ ) / 4f(-)^x [/(+); /(-)]exp -e A (*-*+>) _e A ^-A-)) 



+00 



+00 



/ d/ ( _ ) P x [/ (+) ;/ ( _ ) ]^(/ (+) -/ 1 )K/(-)-/ 2 ) 



(19) 



which coincides with the definition, cq. (|T4|) . 

Further calculations of the two-point distribution function V x (fi, fs) to a large extent repeats the procedure de- 
scribed in detail in the previous paper [l3| for the one-point free energy distribution function. Using the definitions, 
eas. pIH) - (ITT1) . the distribution function, eq.tfTg]). can be represented as follows: 



V x {h,h) = lim E — riDi exp(AL/! + Ai?/ 2 ) / dx 1 ...dx L dy 1 ...dy R &(x 1 ,...,x L ,y R ,...,y 1 ",t) (20) 

A— too £ — * L/lrCl I ^ ] , r 

where 



L+R 



N 



^>(xi,...,x N ;t) = Z(xi) Z(x 2 ) — Z(x N ) = Y[ 
with the replica Hamiltonian 



a=l 



<Pa(0)=0 



exp(-0H N [<h,<h,-,4>N]) (21) 



H 



1 /"* / N N \ 

- / dr (^2[dMr)] 2 ~f3u^2s[Mr)-Mr)] 

2J ° Va=l " a^b J 



(22) 



The propagator ^(x;^, eq. (j2"Tj) . describes N trajectories 4>o,{t) all starting at zero (</> o (0) = 0), and coming to N 
different points {x\, ...,xn} at r = t. One can easily show that ^(xji) can be obtained as the solution of the the 
imaginary-time Schrodinger equation 



with the initial condition 



Here the Hamiltonian is 



*(x;0)=nf =1 J(<) 

l*L 1 N 

H = -5 E 9 ^ ~ 2 K /Z S ( Xa ~ Xb "> 

a—1 a^b 



(23) 
(24) 

(25) 



and the interaction parameter k = /3 3 u. This Hamiltonian describes N bose-particles interacting via the attractive 
two-body potential —n6(x). 

A generic eigenstate of such system is characterized by N momenta {q a } (a = 1, N) which are splitted into M 
(1 < M < N) "clusters" described by continuous real momenta q a (a — 1, ...,M) and having n a discrete imaginary 
"components" (for details see jig. l20l - [24l |): 



la 



q Q - y(«Q + 1 - 2r) 



(r = 1, ...,n a ) 



with the global constraint 



M 



E n « = N 



(26) 



(27) 
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A generic solution ^(x, t) of the Schrodinger equation (|23|) with the initial conditions, Eq. (f2~4]) . can be represented in 
the form of the linear combination of the eigenfunctions \l/q M ^(x): 



N 1 r r 



M=l 



where we have introduced the notation 



|C M (q,n)| 2 *W(x)*W*(0) exp{-^ M (q)t} 



„ MT „ + 00 , oo "I M 



(28) 



(29) 



and 8{k,m) is the Kronecker symbol; note that the presence of this Kronecker symbol in the above equation allows 
to extend the summations over rt a 's to infinity. Here (non- normalized) eigenfunctions are [ToL [2~3| 



N 



^ m) (x)=e n 



V a<b 



1 + in 



sgn(x a - x b ) 
1v a - qv b 



exp 



(30) 



where the summation goes over N\ permutations V of N momenta q a , eg. (I2l)l) . over iV particles x a ; the normalization 
factor 



|CV(q, n) 



Ml j K / \ 1 2 

-r-r |gq - qp - -2{n a - np)\ 



(31) 



and the eigenvalues: 



M 



fa) = 4E?* = ^E ^-^X>i-««) 



2/3 

M i 2 

Er 2 ^3 



24/3 



(32) 



The last term in the above expression provides just the trivial contribution to the selfaveraging part of the free energy 
(discussed in the Introduction) and therefore it will be dropped out of the further calculations. 
Using the definition, eq. (|30l) . one can easily prove that 



tf( M) (0) =N\ 



(33) 



In this way the problem of the calculation of the probability distribution function, cq. (|20[) . reduces to the summation 
over all the spectrum of the eigenstates of the iV-particle bosonic problem, which is parametrized by the set of both 
continuous, {<?i, ?m}j and discrete {ni, Um}', (M — 1, N); (N = 1, oo) degrees of freedom. 



IV. TWO-POINT FREE ENERGY DISTRIBUTION FUNCTION 



Substituting eqs. d28|)- d33J) into ea.(|20|). we get: 

OO 

V x {fi,h) = 1+ Km V { -!)^e XLh+XRh 



(34) 



L+R>1 



L+R M 
M=l a=l Ln a = l 



oo r +oo 



E 



dq a 
2-KK,n a 



M 

8(j2n a , L + R) |C M (q,n)| 51 J Ll fl(cbn) 

a=l 



where 



M 



\Qa — Qp- % -j{n a - Tip) 



|(7 M (q,n)| 2 = J! 

a</3 \Qa ~qp- —{n a +np)\ 



(35) 
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and 



= e ee nn 

■p(L.R) -p(L) -p(R) a=l C=l 



q p (L) - q v ( R ) - ik 
1 v w ~ Q-pW 



n 

a<b 



q v {L) - q r (h) - ik 
q v {L) — q v (L) 



n 



q~(n) - <7„<h) + ik 



q v {R) — q v (m 



dxi...dx L exp i^(g y (t) - ie)x a 

qo<xi<....<xl<x a—1 



x<yR<---<yi<+co 



dy R ...dyi exp[i ^{q v w + if)% 



(36) 



Here the summation over all permutations V of (L + R) momenta {q\, qL+R.} over L "left" particles {x\, xl\ and 
i? "right" particles {ya, are divided into three parts: the permutations of L momenta (taken at random 

out of the total list {qi, qL+n}) over L "left" particles, the permutations of the remaining R momenta over R 
"right" particles, and finally the permutations P^'^ (or the exchange) of the momenta between the group "L" and 
the group "R". Note also that the integrations both over cc a 's and over y c 's in eq. (|36[) require proper regularization at 
— oo and +00 correspondingly. This is done in the standard way by introducing a supplementary parameter e which 
will be set to zero in final results. The result of the integrations can be represented as follows: 



M L R 

7i,fl(q,n) = r( i+R ) exp{ixJ2 n aQa} E 1111 

a=l -p(L,R) a=l c=l . 



q v (L) - q v (R) - ik 



q v {L) — q v (R) 



E 



1 



■p(L) 1^(1) + 1^,(1)) ■■■ (<1~(L) + ••■ + iLiL)) a 



<-) 



(-) 



,(-) 



n 



X! ( + ) / ( + ) 11 IX / ( - I i { ! 1 1 

-p(R) Qyp(R) \Q^~f(R) T 1 V (R)) ■■■\ C 1 V (R) T ■•■ + ^V^' C<d 



l<b 



n 



,(-> 



-,(+) 



,(+) 



7+) _*(+> 



where 

<5 , i ±) = 9q ± «e 

and where we have used the fact that for any permutation of the momenta, eq. (|26[) . one has: 

L+R M 

E qv " = ^2 na<ia 

a—1 a — 1 

Using the Bethe ansatz combinatorial identity [l2j |. 



(37) 



(38) 



(39) 



E 



1 



JY 

QpAipi +<7 P 2 )•••(%>! + ••• + %>«) JJ bL 



1 



JY 



n .\ n 

lla=l'? a a<E> L 



(where the summation goes over all permutations P oi N momenta {gi, ...,qjv}) we get: 



ga + qb + in 
q a + qb 



A I 



L R 



7 L , fl (q,n) = i ( L+B > expj^E^^} E ill! 

q=1 -p(L,R) a =l c=l 



q v {L) - q p (R) - ik 
q vi L) ~ q v < a a) 



1 



n 



Ilo=l 9^,(1) a<b . "+ '/,-, ' 



(-) 



(-i) J 



r„(+) 



w R „(+) AA a ^> +„^> 
1 l c =i e<dl_ ^ u (H) 



n 



(+) 



(R) 



IK 



( + ) 



,(+) 



(40) 



(41) 



Further simplification comes from the following important property of the Bethe ansatz wave function, eg. (1301) . It 
has such structure that for ordered particles positions (e.g. x\ < X2 < ... < xjv) in the summation over permutations 
the momenta q a belonging to the same cluster also remain ordered. In other words, if we consider the momenta, 
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eq. (l26l) . of a cluster a, {qf , q% , q" a }, belonging correspondingly to the particles < Xi 2 < ... < %i na }, the 
permutation of any two momenta q" and q" of this ordered set gives zero contribution. Thus, in order to perform 
the summation over the permutations V^ L ' R ' in eq. tHTj) it is sufficient to split the momenta of each cluster into two 
parts: {qf, Qm \\lm 9n }> where m a = 0, 1, n a and where the momenta qf, belong to the particles 

of the sector " L" , while the momenta q^ +1 ...,q% belong to the particles of the sector " R" . 
Let us introduce the numbering of the momenta of the sector "i?" in the reversed order: 

C -> 9 1 

a * a 

1n a -l -> Q 2 



(42) 



where m a + s a = n a and (s.f. eq. ([26 



q*r = la + y (n a + 1 - 2r) = q a + —{m a + s a + 1 - 2r) 
By definition, the integer parameters {m Q } and {s a } fulfill the global constrains 

M 



a=l 
M 

^ Sq = R 



(43) 



(44) 



(45) 



In this way the summation over permutations 'P^'^ in eg. (1331) is changed by the summations over the integer 
parameters {m a } and {s a }: 



E (• 

■p(L,R) 



M 

n 



^ S(m a + s a , n c 



M M 
S (E m ° ' i ) ^lE Sq ' ^ 



(46) 



which allows to lift the summations over L, R, and {n a } in eq. (|34j) . Straightforward but slightly painful calculations 
result in the following expression (see Appendix): 



f 00 (—'i) M M 



M=l 



00 „ 

E (-i) m " +s "- x 

1 „ i - 



x exp< 



\ m a +s Q — 1 

/ 

Q=lLm Q +s Q >l " v 

24/3 



^ , G(q a ,m a ,s a ) 
dq a ~ — — ; - x 



2-Kn(m a + s a ) 



|(7M(q,m + s)| 2 G M (q, m,s) 



(47) 



where 



and 



|C M (q,m + s)| = |M J_ 

Q </3 |<7q - 9/? - "f {fn a + s a +mp + sp)\ 



g(q a ,m a ,s a ) = 



'(s„ + %qJ- ] ) r(m Q - f g Q (+)) r(l + m Q + s a ) 



2(™«+ s «)r(m Q + s Q + f r(m Q + s Q - f g Q (+)) T(l + m a )T(l + s a ) 

The explicit expression for the factor Ga/ (q, m,s) is given in the Appendix, eq. (|A.17[) . 



(48) 



(49) 



Redefining 



and 



(50) 



2A 2 



(51) 



with 

the normalization factor \Cm (q, m + s)| 2 , eq. ([4*5]) . can be represented as follows: 
|C M (q,m + s)| 2 = l A K + 5 «)- A K + ^)-^ + ^| 2 _ 



a< p \X(m a + s a ) + X(mp + sp) - ip a + ipp 

M 

I [ [2\(m a + a a )] x det 



X(m a + a a ) - ip a + X{mp + sp) + ipp 



(52) 



(53) 



a,/3=l,...,M 



where we have used the Cauchy double alternant identity 



Ua<a( a ct-ap)(b a -bp) 



= (-lj^-^det 



a a - bp 



= 1,...M 



(54) 



with a a = p a — iX(m a + s a ) and b a = p a + iXymp + ap). 

After rescaling, eqs. (f50|) - (|52|) . for the exponential factor in eq.(j47jl we find 



t K 
_ 2^( m « + s a)ql + ^g( TO « + s «) 3 + ^ m ah + Xs a f 2 + ix(m a + s a )q a 

-> -X(m a + s a )p a + ^A 3 (m Q + s a ) 3 + Xm a fx + Xs a f 2 + iXx(m a + s a )p a 
The cubic exponential term can be linearized using the Airy function relation 

rl 



exp 



-X d (m a + s a y 



dy a M(y a ) exp A(m Q + s a) Va 



(55) 



(56) 



Substituting eqs. ([53| . ([55|) and ([56]) into eq. (|47|) , and redefining y a — > y a + p a — ixp a , we get 



V x (h,f 2 ) = lim 1 +J2 



n 



M! 

M=l a=l 



dy a dp a 
2tt 



Ai(y a + p 2 a - ixp a ) (-1) 



m a -\-s a — 1 



m a +s a > 1 



x exp 



{Xm a (y a + fi) + Xs a (y a + f 2 )j G\J-^-, m a , s a 



x detk[(Xm a , Xs ai p a );(Xmp, Xsp, pp)] a p=1M G M (Jj, m, s) 



(57) 



where 



K[[Xm, Xs,p);{Xm', As', p 1 )} 



1 



Am + As — ip + Am' + As' + ip' 



(58) 
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The crucial point of the further calculations is the procedure of taking the thermodynamic limit A^oo. In this limit 
the summations over {m a } and {s a } are performed according to the following algorithm. Let us consider the example 
of the sum of a general type: 



A I 



R(y,p) = lim TT 



(-1)"" 1 exp{Xn a y a } 



P. An, n 



(59) 



where $ is a function which depend on the factors Xn a , p a /X as well as on the parameters n a and p a (which do not 
contain A). The summations in the above example can be represented in terms of the integrals in the complex plane: 



A I 



R(y,p) = lim J] 



dZn 



2i J c sin(7T2 Q ) 



cxp{Xz a y a } 



$( -r, P, An, n 



(60) 



where the integration goes over the contour C shown in Fig. 1(a). Shifting the contour to the position C shown in 
Fig. 1(b) (assuming that there is no contribution from infinity), and redefining z — > z/X, in the limit A — > oo we get: 



A I 



i r ja^p^^j lim ^/'p p z f,\ 

Zm J C i Z a A-s-oo \A A J 



(61) 



where the parameters y a , p a and z a remain finite in the limit A — > oo. 





C 


C 





1 2 3 4 5 

(a) 





) 

1 2 3 4 5 

(b) 



FIG. 1: The contours of integration in the complex plane used for summing the series: (a) the original contour C; (b) the 
deformed contour C'\ 



To perform the summations over m a and s a in eq. (|57|) it is convenient to represent it in the following way: 



(-1) 



A I 



A I 



v x (h,f 2 ) = i + Yl u , 

M=l ' a=l 



n 



+oo 



dy a dp a 
2tt 



Ai(y a +Pa~ ixpa) 



where 



A I 



Sm{p,y; /i,/2) = lim Yl 



a=l 



oo 

{-l) ma+Sa - 1 eTq>{\m a (y a + fx)+\s a (y a + f 2 )} 



.m Q +s Q >l 



M 

n 



r fPa 

vX' mcn Sa 



det K[(Xm a , Xs a ,p a ); (Xmp, \sp,pp)] G M y-r, m, sj (63) 



Sm(p,y; /i,/2) = n 

lim \ TT 

A->oo I A A 



dZ\ a dZ2 a ( 27T« 



o(z 2a ) H o(zi Q ) 



(62) 



The summations over m a and s a in the above expression can be represented as follows 

oo oo oo oo oo 

E (-l)"^- 1 = £ (-ir*- 1 S(s a ,0)+ ^(-l)-- 1 ^^)- £ (-1)^-1 VJ (-l)--i (64) 

m n +s n >l m a - 1 s a — 1 m Q - 1 Scv — 1 

Thus in the integral representation, eqs. (|5"9")) - (j()Tj) . for the function Su(p, y; /i 1/2)5 eq. (j6"3")l . we get 

M 

exp<| 



■-'2, 



>{^la (z/a + /l) + Z2„ (y Q + /b) } 



^(x' "X' ~Y) GM (f' X if) I det ^[( Zl "' 22 "'- Pct ) ; ( z l/3^2/3,P/3)] (65) 
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Taking into account the Gamma function properties, r(z)|i z i_>.o = 1/z and T(l + z)\ 
eq.(|39"]>. and G, eq. (|A.17|) , we obtain 

li m 0^ £lH Z2 « ^ _ (zi a + z 2ct + (z la + z 2et - 

a->oo V A ' A ' A 



(z 2q +ipi ] ){zi a 



and 



lim G 

A— >oo 



/P fl f2 

VA' A ' A 



2 |^o = 1) f° r the factors G, 
(66) 

(67) 



Thus, in the limit A —¥ oo the expression for the probability distribution function, eg. (1621) . takes the form of the 
Frcdholm determinant 



v x (A,h) = i+E Srjr II 



A/ 



M! 

M=l a=l 
dzi a dZ2 a ( 2~K% 



dy a dp a 



2n 



Ai(y Q +Pq - izPa) x 



/ c , (27ri) 2 Vzi 
x exp|zi Q (y a + fi) + z 2a (y a + f 2 
= det[i - A] 



/2m 2m 
{ o{z 2a ) H d(Z! 

\Zi n Z2 



1 



^2 Q + iPa ) 



Z-2, 



zi a - ipi +) 



det 



Zla + Z2a - ip a + Zip + Z 20 + ipp 



(a,0) = l,2,...,M 



(68) 



with the kernel 



A[(z u z 2 ,p);(zi', z 2 ',p')] = f p-Ai(y+p 2 -ixp)(—8{z 2 ) + —5{z 1 ) 

J-oo 2tt V zx z 2 



Z\Z 2 



X 1 



z\ 



z 2 + ip( "> 

1 



z 2 



Z\ — ip( +s > 



expjzi(y + /i) +z 2 (y + f 2 )} 



Z\ + z 2 — ip + Z\' + z 2 ' + ip' 
In the exponential representation of this determinant we get 

oo 

V x {fi,h) = exp[-^-Tri M 

Af=l 



(69) 



(70) 



where 



A I 



Tri M = fi 

a=l 



x exp 



f [ + °° dy a dp a , 2 

J J ^ Al [Va +P a - IXPa 

dz la dz 2a f2ni 2m x 

n 



C ; (2m) 2 \z la 
[z la (y a + fi) + z 2a (y a + / 2 )| 



x 

2;,-/ . 1 
U 



)(> 



2i< 



22a + ip, 
1 



.(-) 



Zla l Pa 



Zla + z 2 a - iPa + Zl a +1 + z 2 a +l + Wa+1 



(71) 



Here, by dehnition, it is assumed that 



'1M+1 — Z U 



(i = 1, 2) and pa/+i = Pi- Substituting 



1 



Zla + Z 2a ~ iPa + Z la+1 + Z 2a+ i + ip a+ l Jo 

into eq. (l7T1) . we obtain 

M 



duj a exp 



- (zi a + Z 2a - iPa + Z\ a+1 + Z 2a +1 + iPa+l) Ua (72) 



p OO y p ~\~ OO ^ cl 

TrA M = / dwi ...dwM TT / / — Ai(y +p 2 + w Q + u a -i -ixp) exp{ip(oj a - w Q _i)} S(p, y; /i, / 2 ) 



(73) 
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where, by definition, loq = lom, and 

dz\dz 2 /2m 



S(p,y,fi,h) 



iz 2 (2m' „. . 2-ni _, . 1 \ / z-< \ / 

M—5{z 2 ) + —S( Z1 ) (l + 1 ) (l , 

i) z \ zi z 2 Z\Z 2 J V z 2 + ip ( >J\ Zi — ipy + > 



Z2 



j c , {2m") 2 V z\ 
x exp|zi(y + /i) + z 2 (y + f 2 )} 

Simple integrations provide the following result: 

S(p,y, fi, h) = % + h) + 6(y + f 2 ) - % + h)6{y + f 2 ) - 6{y + h)6{y + f 2 ) exp{ip(/ 1 - f 2 ) - 2ey} 

+ -L^6(y + f 2 )--^6(y + f 1 ) 
p + te p — ic 

" ~5{y + h)B{h - f 2 ) \l - exp{i(p + ie)(A - f 2 )} 



(74) 



p + ie 



+ -6(y + h)0(f 2 - A) [1 - exp{i(p - ie){h - f 2 )}] 

p — it 



(75) 



According to eq. (fT7|) in what follows we will be dealing with the sector f 2 > fx only. In this case the above expression 
simplifies to 



S(p,y,h,f 2 )\ h>fi = 



p + ie 



p — icJ 



S(y + h) - S(y + h) exp{ip(/ x - f 2 )} 



p — l€ 

+ e(y + f 2 )- 9(y + ft) exp{ip(/ 1 - f 2 ) - 2ey} 
Note that at edge of the sector f 2 > f\ for f 2 = fx + (in the limit e — » 0) 

S{p,V,fi,f 2 )\ h=fi+0 = 2n5(p)S{y + f 2 ) 



(76) 



(77) 



Thus, according to cqs. flTUl) and ([?6"|). the two-point free energy distribution function V x (fi,f 2 ), eq. (fl"4")l . (in the 
sector f 2 > fi) is given by the Fredholm determinant, eq. (j70|) . with the kernel 



A(uj,u)') = Ai(w+o/-/ 2 ) 



dp 



Ai(w + u)' + p 2 - ipx - f 2 ) - Ai(w + w' + p 2 - ipx - /i) exp{ip(/i - f 2 )} 



-oo 27T 
+ °° dp_ 



i(p — ie) 



3xp{ip(w — cj')} + 



Ai(w + ui' + p 2 — ipx + y) — I dy Ai(w + uj' + p 2 — ipx + y) cxp{ip(fi — f 2 )} 
h ' J-h 



x exp{ip(uj — uj')} 
with u>, uj' > 0. 



(78) 



V. THE ENDPOINT PROBABILITY DISTRIBUTION FUNCTION 

Substituting the above result, eqs. (|7T)|) and (|75| . into eq. (fTT|) for the endpoint distribution function one obtains the 
following expression: 



r+°° r f+°° / \-i 

W(x) = / dfF^-f) / dujdw'(l-B_A * u , w (/,x) 

J-oo J JO V 



(79) 
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where 



oo 

Fxi-f) = det[l - B- f ] = exp[- £ - TrB^. 



M=l 

is the GOE Tracy- Widom distribution with the kernel 

B- f (uj,uj') = Ai(u;W - /), (w,u/ > 0) 

and 

'■+ 00 dp 1 



(80) 



(81) 



^LJU)' (/*; ^) 



27T p — it 



Ai'(w + lo 1 + p 2 — ipx — f) exp{ip(uj — u/)} 



dy 



dp 
2^ 



p Ai(ui + uj' + p 2 — ipx + y) exp{ip(uj — a/) } 



(82) 



Using the standard integral representation of the Airy function one can easily reduce the above function Q^^i (/, x) 
to the following sufficiently simple form: 



1 f +oa 



dy 



1 



1 



V 9w duj 

where the function of two variables ^(uj;x) is the simple generalization of the Airy function 



(83) 



(84) 



Here the integration contour C in the complex plane starts at infinity with the argument — 7r/3 and ends up at infinity 
with the argument +7r/3. 

Thus, eas.([79|). (|83l) and (l84l) complete the derivation of the probability distribution function for the directed 
polymer's endpoint. Unfortunately, et present stage the analytical properties of this function are not quite clear. The 
study of this function require the special analysis and it will be done elsewhere. 



VI. CONCLUSIONS 



In this paper the explicit expression for the the probability distribution function of the endpoint of one-dimensional 
directed polymers in random potential is derived in terms of the Bethe ansatz replica technique. The result obtained, 
eqs. (|79[) - (l84l) . looks quite similar to the one derived in terms of completely different method in which the maximal 
point of the Airy 2 process minus a parabola is considered [l4T - [l6| . Unfortunately, at present stage the final expression 
for the probability distribution function obtained both here and in Refs. [Ml-llo] is rather sophisticated so that the 
study of its analytical properties would require special efforts. Hopefully this problem will be solved in the near future. 

One more conclusion of the present study is that the approach used, namely the Bethe ansatz replica technique, 
ones again (following the works [Tol - [l3| ) has demonstrated its efficiency. Hopefully it will also be fruitful for the 
studies of more serious problems in this scope, such as joint statistical properties of the free energy fluctuations at 
different times. 
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Appendix 



In terms of the parameters {m a } and {s a } the product factors in eq. (|41[) are expressed as follows: 



n<?> 



M m D 



(-) 

R Ms 



nn*? ( - ) 



(A.l) 
(A.2) 



n 

a<b 



M m a 

n n 

a— 1 l<r<r' 



«(-) 



•t?2 (-) + IK 



>< n nn 

l<a</3 r=l r'=l 



«(-) 



q? [ ~ ] +q 



0<-> 



(A.3) 



n 



+*(+) 



M s Q 

n n 

q— 1 l<r<r' 



g*? (+) +9*° (+) -« 



n nn 



l<a</3 r=l r'=l L g 



„*«(+) , */3 ( + ) 



ko(+) ,„*0 (+) 



(A.4) 



nn 

a=l c=l 



q v (h) — q v (R) — ik 
q v ( a L> ~ q v ( c R) 



M ( m a sp 

n nn 

l<a</3 Lr=l r' = l 



V * P 

— q r , — in 
q? + q*r' 



nn 



q*r 



6 



g*7 - ql 



M m a s a 

nnn 

a=l r—1 r' — l 



g" — q r , — ik 
q? - q*r> 



(A.5) 



Substituting eqs. (IA.l[) - (|A.5[) into eq. (|4"Tj) . and then substituting the resulting expression into eq. (|3~4")) we obtain eq. (jT7)) 
where 



n 



g*° (+) +g*°, (+) -m 



nn 



g" — g r , ~ ik 
q? ~ q*r> 



(A.6) 



and 



^=nn 



ct(-) i ,8 ( _ ) , ■ ' 
?r + V + ZK 



nn 



? *a(+) +(? *£ 



(+) 



m S£ 

nn 

r— 1 r' — l 



g? - - IK 



r=lr' = l if" /—I/' i '/ -:■'/' ■ /-I.,'- J 9r + 9 

The product factors in eq. (|A.6[) can be easily expressed it terms of the Gamma functions: 



nn 



q r - q 1 ;, - ik 

q*r - <& 



(A.7) 



n^ (_) = fl[fc (_) - ^{m a + Soi + l)+mr\ = (my 



p[ 1 s a —m a ig a 

1 2 2 



1 I 2 2 



(A.8) 



i[q*r {+) = n[^ (+)+ Y( m « +s « +i )-H = 



rfi- 



(A.9) 



n 



qf- ] + g?, (-) + 
q^ ] + q?S~ ) 



= 2~ (r ' 



n 



(-) ' 



2 -(««.-i) 



r( Sc 



r 1 



2iq Q 



r l- 



r( 1 — m Q + 



2i<;,, 



nn 

r=l r' = l 



— <?*" — IK 
9? - q*r' 



r(i + m a + s a ) 

r(i + m a )r(i + s Q ) 



Substituting the above expressions into eq. (|A.6[) and using the standard relations for the Gamma functions, 

r(*or(i-«) 



sin(7rz) 
r(l + z) = zT(z) 



v^T(i + 2z) 



for the factor Q, ea. (|A.6|) . we get 

G(q a ,m a ,s a ) 



2 2z r(i + z) 
r( Sa + ^ Q ( - ) )r(m Q -^ Q (+))r(i 



2(«<.+»«.)r(m Q + s a + %qj-)) T(m a + s a ~ f 9a (+)) r(l + m a )r(l + s a ) 
Similar calculations for the factor <7a/3, eq. (|A.7|) yield the following expression 



£ Q/3 (q, m, s) 



1+ TOa+m ,- Sa - s . _i (ga (- )+ ^(-) 



m a -\-mp — s a —S0 



r 



S-) + n a (-)\ 



m a —mi3—s a —Sfi 



( + ) 



771^ +^i/3+-5 a — Sfl 



+ j(?aW+3fl( + )) 



9/3) 



2 



-m a +mg+s a — 



-9/3) 



1 + 



-m a +mp+s a +sp 



2 










i 


2 




K 


m^+m^ — s Q 




i 


2 




K 



9/3) 



9/3) 



1 T 2 



:(«fa-«|B) 



i 9a 



9/3) 
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